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ABSTRACT 



Several dynamical processes may induce considerable electric currents in the atmospheres of magnetic chemically peculiar (CP) stars. 
The Lorentz force, which results from the interaction between the magnetic field and the induced currents, modifies the atmospheric 
structure and induces characteristic rotational variability of the hydrogen Balmer lines. To study this phenomena we have initiated a 
systematic spectroscopic survey of the Balmer lines variation in magnetic CP stars. In this paper we continue presentation of results 
of the program focusing on the high-resolution spectral observations of AOp star 9 Aur (HD40312). We have detected a significant 
variability of the Ha, HyS, and Hy spectral lines during full rotation cycle of the star. This variability is interpreted in the framework 
of the model atmosphere analysis, which accounts for the Lorentz force effects. Both the inward and outward directed Lorentz forces 
are considered under the assumption of the axisymmetric dipole or dipole+quadrupole magnetic field configurations. We demonstrate 
that only the model with the outward directed Lorentz force in the dipole+quadrupole configuration is able to reproduce the observed 
hydrogen line variation. These results present new strong evidences for the presence of non-zero global electric currents in the 
atmosphere of an early-type magnetic star. 
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1. Introduction 

The atmospheres of magnetic chemically peculiar (CP) stars dis- 
play the presence of global magnetic fields ranging in strength 
from a few hundred G up to several tens of kG (lLandstreetl 
|2001|) . In the contrast to the complex, localized and unstable 
magnetic fields of the cool stars with an external convective 
envelope, magnetic fi elds of CP stars are organized at a large 
scale, roughly dipola r dLandstreetl2 00 1 ) or low-order multipolar 
(Bagnuloetd][2002) geometries, which are likely stable during 
significant time intervals. In fact, these stars provide a unique 
natural laboratory for the study of secular evolution of global 
cosmic magnetic fields and other dynamical processes which 
may take place in the magnetized plasma. In particular, the slow 
variation of the field due to decay changes the pressure-force bal- 
ance in the atmosphere via the induced Lorentz force, that makes 
it possible to detect it observationally and establish a number of 
important constraints on the plausible scenarios of the magnetic 
field evolution in early- type stars. 

Several physical mechanisms have been suggested as the 
sources of non-zero Lorentz force in the atmo spheres of CP 
stars: the global field distortion and evolution (St epienlll978t 
lLandstreettl987HValvavin et al.l2004l) . the global drift of charge 
atmospheric particles un der the influence of r adiative forces 
dPeterson & Thevsl 1 198 lb . ambipolar diffusion (LeBlan c et al.l 
119941) . Phenomenological atmosphere mo dels w i th the L orentz 
force incl u ded w ere also presented by iMadeil (ll983alTbh and 
Carp enter! d!985l) . In the light of these discussions it becomes 
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clear that the study of the Lorentz force in the atmospheres of 
CP stars is of fundamental importance for understanding the 
nature of intrinsic microscopic processes in magnetized atmo- 
spheric plasma. 

The aforementioned studies suggested that the magnetic 
forces may lead to significant differences between the atmo- 
spheric structures of magnetic and non-magnetic stars. In some 
cases, the Lorentz force may noticeably change the effective 
gravity and influence formation of the pressure-sensitive spectral 
features, especiall y the hydrogen Balmer lines. Some of the H/j 
photometric data dMadei et alj|1984t iMusielok & Madejlll~98"8l) 
can be considered as an evidence for the presence of non-force- 
free magnetic fields. S pectroscopy of the hydrogen lines also 
points in this direction. iRrollI d 19891) found variability with am- 
plitudes more than 1% in the Balmer lines in several magnetic 
stars. He showed that at least part of this variability can be at- 
tributed to the presence of a non-zero Lorentz force in stellar 
atmospheres. 

Unfortunately, apart from the study by Kroll, who presented 
low-resolution spectroscopic observations of the Balmer line 
variability with rotation phase in only a few CP stars, there have 
been no other systematic spectroscopic surveys of the Balmer 
line variability. Such a situation is due to the fact that the ob- 
servational aspect of the problem is fairly complex, and com- 
prehensive understanding can be obtained only with the help of 
high-precision, high-resolution spectroscopic observations. Very 
weak variations in the Balmer line profiles (for the majority 
of stars it is 2% or less) require reconstruction of the strongly 
broadened spectral features with an accuracy of ~ 0. 1 %. Until re- 
cently, such a precision could not be reached by high-resolution 
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spectroscopy. However, the situation has been improved signif- 
icantly with the development of highly stable fibre-fed spec- 
trographs that made it possible to carry out these studies at a 
more advanced instrumental level. Taking this in to account and 
following the pioneering work by iRrollI (11989) we have initi- 
ated a new spectros copical search of hydrogen line variability 
dValvavin et al.ll2005l) . 

In this paper we present the phase-resolved high-resolution 
observations of one of the brightest weak-field {B as 1 kG) mag- 
netic CP star Aur. We detected significant variation of the 
Balmer line profiles and interpreted it in terms of the non-force- 
free magnetic field configuration. We argue that chemically over- 
abundant spots in the atmosphere of 6 Aur can not produce the 
observed variability. 

In the next section we describe observations and methods of 
spectral processing. Variability of the Ha, HB, and Hy lines is 
illustrated in this section. Sect. 3 introduces the model which 
we have employed to interpret the observations. In Sect. 4 we 
perform calculation of our stellar model atmosphere with the 
Lorentz force. Results are presented in Sect. 5 and summarized 
in Sect. 6. General discussion is presented in Sect. 7. 



2. Observations 

Aur (HD40312) is a broad-lined AOp st ar with a relatively 
weak («lkG) dipolar magnetic field (see IWade et all 12000). 
During rotation the star shows equatorial (phases 0.25, 0.75) as 
well as polar regions (phases 0.0, 0.5) of its magnetosphere . As 
follows from the previous studies (e.g. Valvavin et al. 2004]), the 
maximum atmospheric perturbation by the Lorentz force is ex- 
pected to be observed at the equatorial plane and is nearly zero at 
the polar regions. This makes it possible to estimate the magnetic 
force term by analyzing the differences between the Balmer line 
profiles obtained at different rotation phases. 

The observations were carried out with the BOES echelle 
spectrograph installed at the 1.8 m telescope of the Korean 
Astronomy and Space Science Institute. The spectrograph and 
observational procedures are described bv lKim et alJ fcOOO). The 
instrument is a moderate-beam, fibre-fed high-resolution spec- 
trograph which incorporates 3 STU Polymicro fibres of 300, 
200, and 80 /mi core diameter (corresponding spectral resolu- 
tions are A/AA = 30000, 45 000, and 90000 respectively). The 
medium resolution mode was employed in the present study. 
Working wavelength range is from 3500 A to 10000 A. High 
throughput of the spectrograph at 4100-8000 A wavelength 
range and its high stability make it possible to obtain spectra 
of the Balmer lines with an accuracy of about 0.2 - 0.3%. 

Twenty spectra of 6 Aur were recorded in the course of about 
20 observing nights from January 2004 to April 2005. Typical 
exposure times of a few minutes allowed to achieve S /N ~ 150. 
Table [T] gives an overview of our o bservations. Throu ghout this 
study we use ephemeris derived by Wad e et al] (|2000): 

7D = 2450001.881 + Ex 3.61860, (1) 

where the reference time corresponds to the negative ex- 
tremum of the longitudinal field variation. The spectral reduc- 
tion was carried out using the image processing program DECH 
(Galazutdinov 2000) as well as MIDAS packages. The general 
steps are standard and include cosmic ray hits removal, elec- 
tronic bias and scatter light subtraction, extraction of the spectral 
orders, division by the flat-field spectrum, normalization to the 
continuum and wavelength calibration. 



Table 1. Observations of 9 Aur. The first column gives the num- 
ber of observation, the second column lists the Julian date (JD) 
an d the last one is th e rotation phase calculated with ephemeris 
bv lWade et alJd2000h . 



No. 


JD 


Rotation Phase 


1 


2453015.0238 


0.681 


2 


2453015.9949 


0.950 


3 


2453020.0157 


0.061 


4 


2453038.0886 


0.056 


5 


2453039.1761 


0.356 


6 


2453039.9192 


0.562 


7 


2453040.9137 


0.837 


8 


2453042.9287 


0.393 


9 


2453046.0285 


0.250 


10 


2453340.0242 


0.496 


11 


2453341.0364 


0.775 


12 


2453341.3455 


0.861 


13 


2453343.0469 


0.331 


14 


2453354.0668 


0.377 


15 


2453354.3565 


0.456 


16 


2453356.0911 


0.936 


17 


2453356.3045 


0.995 


18 


2453457.9726 


0.091 


19 


2453458.0466 


0.111 


20 


2453458.0834 


0.121 




Fig. 1. Examples (from top to bottom) of the initial image of 
echelle orders, its 2-D continuum and normalized image. 



Normalization to the continuum deserves some additional 
comments. In order to obtain uniformly reconstructed contin- 
uum in all the spectra we applied the following technique. After 
the flat-fielding procedure all extracted spectral orders of in- 
dividual frames were merged into new 2-D (4000x75) images 
where the echelle orders (75 orders in total) where consequently 
placed as image rows. Then, with the aid of the median and 
Gaussian filters, we identified and cut out all narrow spectral 
features in all spectral orders. Finally, ignoring the Balmer line 
regions we fitted all the images using 2-D cubic spline function 
and created 2-D continuum images. Normalized spectra were 
produced by the division of the initial images by their corre- 
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Fig. 2. Mean spectrum (the upper plot) and standard deviation fl- 
at the H/3 and Hy profile of Vega spectrum. 



sponding 2-D continuum images. Examples of the initial image, 
2-D continuum derived from it and normalization results are pre- 
sented in Fig.Q] 

Our analysis showed that such a technique makes it possi- 
ble to achieve a required high accuracy and stability of the con- 
tinuum reconstruction in those places of the observed spectra 
where the spectrograph's response function can be described by 
a monotonic low-order polynomial 2-D function. For the BOES 
the most appropriate region is from 4150 to about 5100 A that 
makes our analysis reliable at the Hy and H/3 Balmer lines. 
Accuracy of the continuum normalization around these lines 
is estimated to be approximately 0.1-0.2%. In order to illus- 
trate this conclusion, here we present results of test observa- 
tions of Vega which were carried out in different nights in 
the period from earl y March to May of 2006. It is established 
dPeterson et al.ll2006b that this standard star has the inclination 
of the rotational axis i = 4.5° ± 0.33°. Such a small inclina- 
tion minimizes the probability to detect any physical variability 
in the Balmer profiles of Vega's spectrum allowing to use this 
star as a standard in our study. Results of the tests are presented 
in Fig. |2] where the upper plot illustrates mean spectrum of the 
star obtained by averaging the observed spectra. The lower plot 
is the dispersion of the result obtained as standard deviation of 
the individual spectra from the mean. As can be seen from the 
behavior of the dispersion function, all instabilities of the con- 
tinuum reconstruction lie typically below 0.15 - 0.2%. 

Also, we would like to note that during these observations 
the spectrograph's configuration has been changed many times 
by service requirements. Nevertheless, the response function of 
the spectrograph showed very good stability despite the fact of 
these re-configurations. Analyzing observations of another pro- 
gram Ap/Bp stars which do not show any significant variabil- 
ity of the Balmer profiles we also concluded that the response 
function does not change its shape for the last three years (2003- 
2006). 

For different reasons, the stability around the other Balmer 
lines is not so good compared to the Hy and H/3 regions. For 
example, the He region is characterized by the presence of a 
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Fig. 3. The standard deviation cr at the Ha, H/3 and Hy. The pho- 
ton noise background of cr as 0.3% is subtracted. 



narrow slope of the response function that makes it difficult to 
reconstruct continuum around this line with the necessary accu- 
racy. However, in this study we decided to include the Ha line 
also as an illustration. 

Taking all the above conclusions into account we searched 
the Balmer line profile variations using normalized spectra at the 
Ha, H/3, and Hy regions. For each of the spectral intervals we 
examined the standard deviation <x from the mean in the profiles 
during full rotation cycle of the star. The standard deviations as 
function of wavelength are presented in Fig. [3] Analyzing spec- 
tral regions with low lines density we found that the standard 
deviation due to the photon noise and inaccuracies in the spectra 
processing lie below the level of 0.3%. Any deviation above this 
value indicates the presence of significant intrinsic variability of 
Balmer line profiles. 

The standard deviations at the H/3 and Hy reve al the char- 
acteristic fingerprint first described by iRrollI (Qj89) as the im- 
pact of the Lorentz force: the amplitude rises in the wings and 
drops to the line center. This picture is less clear for the Ha 
line due to the above reasons and strong distortion of the pro- 
file by the telluric molecular absorptions. (Here we note that the 
spectra obtained at different dates on the time base of about one 
year are strongly Doppler-shifted relative to each other due to 
the Earth's movement. In this connection our study requires re- 
binding the spectra to the rest wavelengths and, as a result, all 
the telluric lines become randomly shifted in the individual re- 
bined spectra that strongly complicates analysis of the Ha pro- 
file). Nevertheless, the effect is seen in the red wing of the line 
(see Fig. [3}- 

The narrow features in the standard deviations (H/3 and Hy 
lines, Fig. [3} come from metal lines located in wings of the 
hydrogen line profiles. Strong variability of these lines indi- 
cates the presence of an inhomogeneous surface distribution of 
corresponding chemical abundances. This distorting factor con- 
tributes to an additional noise but it can be isolated and separated 
from the broad spectral changes in the hydrogen lines. 
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3. Model 

3.1. General equations and approximations 

In order to model the found variation of the Balmer line profiles 
we assume that this variability is caused by the Loren tz force 
and follow approaches outlined bv lValvavin et al.l d2004l) . 

In the presence of the Lorentz force term the hydrostatic 
equation reads: 



VP, 
and 



total 



: Pg + fh, 



(2) 



(3) 



/l = -JXB, 

c 

where P t otai = Pgns + Prad denotes the total pressure, P gas and P lac | 
are the gas and radiation field pressures respectively, p is the gas 
density, g is the surface gravity, j and B are the surface currents 
and magnetic field vector respectively, and c is the speed of light. 
The electric currents are determined using the Ohm's law: 



j = AE\\ + A ± E ± + A H 



BxE ± 

B ' 



(4) 



where E\\ and E± are electric field components directed along 
and across magnetic field lines respectively, A is the electric con- 
ductivity in the absence of the magnetic field, A ± is the electric 
conductivity across the magnetic field lines and An is the Hall's 
conductivity. 

To simplify solution of the problem, we consider poloidal 
surface magnetic field geometry. The non-force-free term giv- 
ing rise in the Lorentz force is described via the induced electric 
field [e.m.f.). This configuration can be justified in the context 
of the magnetic field evolution. For example, distortion of ini- 
tially force-free configuration of the magnetic field and respec- 
tive distribution of the electric field over the stellar surface may 
be produced by the global field decay/generation or by any other, 
more complicated scenario related, for instance, to generation of 
internal toroidal fields by differential rotation etc. In the present 
study we do not consider these details and generalize the prob- 
lem making no assumption about the origin of the Lorentz force. 
Following the approach developed in previous studies and tak- 
ing into account the fact that the surface magnetic field of the 
majority of CP stars may approximately be described by low- 
order axisymmetric poloidal fields, we restrict our model con- 
sidering only azimuthal geometry of the induced electric field. 
In this case the equations introduced above can be implemented 
in a 1-D stellar model atmosphere code. 

Finally, the complete set of assumptions used in our mod- 
elling can be summarized as follows: 

1. The stellar surface magnetic field is axisymmetric and is 
dominated by dipolar or dipole+quadrupolar component in 
all atmospheric layers. 

2. The induced e.m.f. has only an azimu thal component, sim- 
ilar to that described bv lWrubell (119521) . who considered de- 
cay of the global stellar magnetic field. In this case the dis- 
tribution of the surface electric currents can be expressed by 
the Legendre polynomials P l n {n), where n — 1 for dipole, 
n — 2 for quadrupole, etc., and // = cos 9 is the cosine of the 
co-latitude angle 9 which is counted in the coordinate system 
connected to the symmetry axis of the magnetic field. 

3. The atmospheric layers are assumed to be in static equilib- 
rium and no horizontal motions are present. 

4. Stellar rotation, Hall's currents, ambipolar diffusion and 
other dynamical processes are neglected. 



Taking into account these approximations and substituting 
Eq. (0 and Eq. © into Eq. (|2]), we can write the hydrostatic 
equation in the following form 



dftotal 

dr 



= -pg± 1 -A ± Y J c n P 1 MZ B{ e ) - 



(5) 



Obtaining this equation we used the superposition principle for 
field vectors and the solution of Maxw ell equa t ions f or each of 
the multipolar components following Wrubel (1952). We also 
suppose that E ± B. Here c„ represents the effective electric 
field generated by the «-th magnetic field component at the stel- 
lar magnetic equator and Be is the horizontal field component. 
The signs "+" and "— " refer to the outward and inward directed 
Lorentz forces, respectively. Equation © can be solved if the 
plasma conductivity A± is known. 

We note that the values of c„ are free parameters to be found 
by using our model. These values represent the fundamental 
characteristics which can be used for building self-consistent 
models of the global stellar magnetic field geometry and its evo- 
lution. Thus, an indirect measurement of these parameters via 
the study of the Lorentz force is of fundamental importance for 
understanding the stellar magnetism. 

3.2. Calculation of plasma conductivity 

Calculation of the electric conductivity Ax can be carried out 
using the Lorentz collision model where only binary collisions 
between particles are allowed. The detaile d description and basi c 
relationships of this approach is given in IValvavin et all d2004l) . 
Improving the latter work, here we calculate the electric conduc- 
tivity including all available charged particles. The precise cal- 
culations require direct evaluation of the effective stopping force 
acting between two charged particles of types k (background par- 
ticle) and i (test particle): 



1 eff ~ 



47r(Z,e) 2 (Z it e)VlnA 



fik 



(6) 



Here and Vjt,T are the concentration and the most probable 
thermal velocity of particles of type k, Z, and are charges of 
the particles, M is the reduced mass and In A is the classical 
Coulomb logari thm. The Chandrasek har function fix) has the 
following form (IChandrasekharl 1 942b : 



/(*) = — f* 



Jo 



dx — xe 



(7) 



In the case of the interaction between charged and neutral parti- 
cles we use the elastic collision cross-sections, 



2 m k 



-2/3 



(8) 



where «o is the Bohr radius and m p is the proton mass. 

Generally, at the uppermost atmospheric layers where the 
cyclotron frequency of the conducting particles is much higher 
than their mean free-path times, effects of magnetoresistivity on 
A ± redire ct induced horizonta l electric currents into Hall's cur- 
rents (see Val yavin et al.ll2004l) which are ignored in our model. 
In the deepest layers, the effects of the magnetic field influence 
on the conductivity become negligible and electric currents can 
be described by ordinary Ohm's law for the non-magnetic case. 
In the intermediate atmospheric layers, however, the situation 
becomes more complicated. Under the influence of horizontal 
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magnetic field, the contribution of ionized particles to the con- 
ductivity is strongly stratified by the magnetic fiel d that produces 
the b imodal shape of the Lorentz force term dValvavin et al.l 
12001 . 

Computations of the electric conductivity under the afore- 
mentioned assumptions enables us to present the solution via un- 
known induced e.m.f. (c„) and magnetic field d etermined from 
obser vations. A preliminary analysis (see also IValvavin et al.l 
2004) showed that amplitudes of the pressure-temperature mod- 
ification with the Lorentz force depend mainly on the induced 
e.m.f. Magnetic field strength also influences the amplitudes 
but mainly affects the positions of the perturbed bimodal area 
on the Rosseland optical depth scale. For the very small mag- 
netic fields the perturbation of the atmosphere is limited to small 
Rosseland optical depths above the region where the wings of the 
Balmer lines form. Increasing the magnetic field strength shifts 
the perturbation toward the deeper regions, which contribute to 
the formation of hydrogen lines. In this process the two maxima 
of the bimodal atmospheric perturbations pass one after another 
through the zone of Balmer line formation, giving rise to max- 
imum perturbations to the Balmer line formation around hori- 
zontal local magnetic field B e of 400-500 G and 10000 G with a 
significant gap between Be — 1 kG and 6 kG. In terms of global, 
nearly dipolar magnetic fields these areas of the most effective 
contribution of the horizontal magnetic field to the Balmer line 
formation correspond to approximately 1 kG and 20 kG mag- 
netic stars. Taking into account that the majority of magnetic 
stars have magnetic fields weaker than 10 kG we predict that 
the maximum amplitudes of the Balmer line variation due to the 
Lorentz force are expected to be among magnetic Ap/Bp stars 
with magnetic fields between 0.4 kG and 1-2 kG. The Aur is 
one of such stars. 



3.3. Model atmospheres with Lorentz force 

Our calculations were carried out wit h the stellar model atmo- 
sphere code LLmodels developed bv lShulvak et alJ f2004). The 
code is bas ed on the mod ified ATLAS9 dKuruczl ll993a) and 
ATLAS 12 dKuruczlll993bl) subrout ines as well as o n the spec- 
trum synthesis package described by Tsymbal ( 1996). LLmodels 
is written in Fortran 90 and uses the following general approx- 
imations (which are typical for many 1-D model atmosphere 
tools): 

- the plane-parallel geometry is assumed; 

- the Local Thermodynamic Equilibrium (LTE) is used to cal- 
culate the atomic level populations for all chemical species; 

- the stellar atmosphere is assumed to be in a hydrostatic equi- 
librium; 

- the radiative equilibrium condition is fulfilled. 

The code incorporates the so-called line- by-line (LL) metho d 
of the bound-bound opacity calculations ( Shul yak et alj |2004). 
This technique allows us to account for individual stellar abun- 
dance pattern, which is important in the present study of 9 Aur 
because this star shows non-solar abundances and an inhomoge- 
neous horizontal distribution of some chemical elements. In this 
case accurate treatment of lines opacity in model atmospheres 
is needed to ensure correct calculation of model structures with 
individual abundance patterns. 

The new module was written to compute electric conductiv- 
ity and magnetic pressure. Since the effective gravity is a func- 
tion of atmospheric depth, additional changes in the solution of 
the hydrostatic equation and in the mass correction routines were 
made. The hydrostatic equation is solved and presented in terms 




4000 



5000 6000 
Wavelength (A) 



7000 



8000 



Fig. 4. Comparison of the observed (symbols) and computed 
(solid line) energy distribution of 6 Aur. Theoretical model cor- 
responds to r e ff = 10 400 K, log g = 3.6, and has been convolved 
with a FWHM = 20 A Gaussian filter. 



of monochromatic optical depth scale 75000 as an independent 
variable. At each iteration the code calculates electric conduc- 
tivity in all atmospheric layers using all available charged and 
neutral plasma particles. The conductivity is then used to evalu- 
ate magnetic contribution to the effective gravity and to execute 
mass correction procedure. 

As can be seen from Eq. ((5]), in the case of the outward di- 
rected Lorentz force, there is some critical value of c„ which may 
produce unstable solution. Such models can not be considered in 
the hydrostatic equilibrium approximation introduced above and 
were assumed to be non-physical in our calculations. Thus, for 
each set of models, we adopted c„ values as to ensure static equi- 
librium. 

In addition, the following calculation settings have been 
used: the atmosphere is sliced into 72 layers equally 
spaced along the Rosseland optical scale height tr oss , from 
tross = -6.875 to tr oss = 2.0. The number of frequency points 
used for the flux integration procedure was 495 000 in the 
500-50 000 A spectral region. The initial atomi c line list taken 
from VALD (Vienna At omic Line Database) ( Pis kunov et al.l 
119951; lKupkaetal.l [l999) contains information about 21.6 mil- 
lion atomic lines, including lines originating from the predicted 
energy levels. This line list was used as an input for the preselec- 
tion procedure in the LLmodels code. We adopted the selection 
threshold { v /a y > 1%, where a v and { v are the continuum and 
line absorption coefficients at the given frequency v. The pres- 
election enabled us to reduce the total number of lines used for 
the line opacity calculations to about 525 000. 

4. Results 

4.1. Model atmosphere parameters of 6 Aur 

The model atmosphere parameters, \ogg and T e g, were 
determined using s pectrophotometric energy distribution 
dAdelman et all 1 1987b and theoretical fit of the WB and Hy 
line profiles. For this purpose we have chosen observations at 
phase 0.995 (polar region) where the influence of the Lorentz 
force is small and the stellar atmosphere is expected to be 
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Fig. 5. Comparison of the observed and computed HB (upper 
panel) and Hy (lower panel) profiles. Thick line shows observa- 
tions at phase 0.995, thin line corresponds to synthetic spectra. 

Table 2. The surface-averaged abu ndances (log jV /^Vtotai) of 
9 Aur for the four rotation phases dKuschn ig 1998) and aver- 
aged over the whole stellar surface (last column). 



Element 


0.00 


0.25 


0.50 


0.75 


average 


He 


-2.32 


-2.32 


-2.40 


-2.32 


-2.34 


Mg 


-5.28 


-5.27 


-5.35 


-5.50 


-5.35 


Si 


-3.35 


-3.27 


-3.09 


-3.22 


-3.23 


Ti 


-7.52 


-7.61 


-7.67 


-7.61 


-7.60 


Cr 


-5.14 


-5.35 


-4.99 


-4.75 


-5.06 


Mn 


-5.50 


-5.64 


-5.49 


-5.34 


-5.49 


Fe 


-3.86 


-3.83 


-3.63 


-3.69 


-3.75 


Sr 


-8.43 


-8.43 


-8.36 


-8.19 


-8.35 



close to the standard non-magnetic model structure. The chem- 
ical composition as well as the p rojected rotation al velocity 
v sin / = 55 kms 1 were taken from lKuschnigl (fi998). who used 
multi-element Doppler Imaging technique to derive surface 
maps for eight elements in 9 Aur. Table [2] gives abundance pat- 
terns for four representative rotation phases. Synthet ic Balmer 
line p rofiles were calculated using the Synth program dPiskunovl 
|1992|) . The program incorporates recent improvements in th e 
treatment of the hydrogen line opacity (Barkl em et al.1 12000). 
The stellar energy distribution and Balmer lines are best approx- 
imated with the following parameters: T e g = 10400 ± 300 K, 
logg = 3.6 ± 0.05. Note, that such a high accuracy of the de- 
termined parameters is just an internal accuracy obtained from 
our technique which we used to fit the data. Real parameters 



may be slightly different from the obtained ones due to various 
systematic error sources, but this does not play a significant 
role in our study. Comparison of the observations and model 
predictions are presented in Fig. [4] and Fig. [5] Fundamental 
parameters of 9 Aur suggest that this star is significantly 
evolved from the ZAM S. From the compar i son w ith theoretical 
evolutionary tracks Kochukhov & Bagnulo (2006) find 9 Aur to 
be at the very end of its main sequence life. 

Recent studies b y iKochukhov et al.l ([2005 ) and 
iKhan & Shulvakl d2006allbh showed that the effects of Zeeman 
splitting and polarized radiative transfer on model atmosphere 
structure and shapes of hydrogen line profiles are less than 0. 1 % 
for magnetic field intensities around 1 kG and thus can be safely 
neglected in the present investigation. 

4.2. Magnetic field geometry 

To calculate the Lorentz force effects it is essential to spec- 
ify the magnetic field geometry (see Eq. (O). The first longi- 
tu dinal magnetic field measu rements were obtained for 9 Aur 
bv lBorra & Landstreetl (fi980) using HB photopolarimetric tech- 
nique. The authors observed a smooth single-wave (Z? z ) variation 
with rotation phase and concluded that it is pr obably caused by a 
dipole inclined to the rotation axis of the star. lWade et alj (12000) 
presented high-precision longitudinal field measurements of this 
star. T hey confirmed and improved results of lBorra & Landstreetl 
(1980). In particular, deviation of the (B z ) variation from the 
p urely sinusoidal m agnetic curve has been found (see Fig. 12 
in lWade et alj|2000l) . 

Following these works, we have approximated magnetic 
field topology of 9 Aur by a combination of the dipole and ax- 
isymmetric quadrupole magnetic components. We have also as- 
sumed that the symmetry axes of the dipole and quadrupole mag- 
netic fields are parallel. Thus, the model parameters include the 
polar strength of the dipolar component Ba, relative contribution 
of the quadrupole field B q /B^, magnetic obliquity B, and inclina- 
tion angle i of the stellar rotation axis with respect to the line of 
sight. The last parameter is best estimated independently, from 
the usual oblique rotator relation connecting stellar radius, rota- 
tion period and v sin i. Hipparcos parallax ti = 18.83 + 0.81 mas 
and r eff = 10400 + 300 K yield R = 5.1 + 0.4 R Q . This leads 
to i = 51° ± 6°, which is in good agr eement with =50°- 
60° derived in Doppler imaging studies dRice & Wehlaul [l991: 
lKuschnidll998l) . 

The remaining free parameters of our model for the magnetic 
field geometry were dete rmined with the least-squares fit of the 
obser ved (B z ) variation dBorra & Landstreetl 1 19801: IWade et al.1 
2000). We have examined "positive" (quadrupolar field has two 
positive and one negative poles) and "negative" (the two neg- 
ative and one positive poles) quadrupolar field configurations. 
We found that the lowest^ = x 2 / v corresponds to the nega- 
tive quadrupolar configuration with B q /Ba ~ -2 (xl = 1.14, 
B& = 1 .4 kG, B = 78°), while the purely dipole model gives 
^ = 2(B d = 1.3kG, J 8 = 73°). 

4.3. Effects of horizontal abundance distribution 

In order to distinguish effects of the magnetic pressure from the 
ones of the abundance distribution, we calculated model atmo- 
spheres for each of the four representative phases using individ- 
ual abundances from Table|2] The synthetic profiles of the hydro- 
gen lines as well as the standard deviation of the profiles due to 
the variable chemical abundances at the listed phases were then 
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Fig. 6. The standard deviation cr at the Ha, WB and Hy lines. 
Observations are shown by solid lines. The dashed lines illustrate 
variability due to the horizontal abundance inhomogeneities. 

calculated. The results are presented in Fig. [6] As one can see, 
effects due to a non-uniform abundance distribution are totally 
different from the observed one: the strongest variability occurs 
in the line core whereas the line wings are not affected much. 
This allows us to conclude that the chemical spots do not lead to 
the observed Balmer line var iations. This fact is in concordance 
with the results obtained by iRrollI d 19891) . who showed that the 
found shape of the Balmer line variability can not be reproduced 
by metallicity or temperature variations, but can be considered 
in the frame of changes in the pressure structure of the stellar 
atmosphere. Finally, the average abundances from Table [2] (last 
column) were used in the rest of model atmosphere calculations. 



4.4. The Lorentz force 

In order to fit the observed variations of the line profiles, both 
the inward and outward directed Lorentz forces were examined 
through the model atmosphere calculations. The actual magnetic 
input parameters of computations with LLmodels code include 
the sign of the Lorentz force, magnetic field modulus B, and the 
product of two sums Y> Cnfl £ Following the procedure de- 
scribed bv lValvavin et al.l (120041) . we take the latter two param- 
eters to be disk-averaged at the individual rotation phases. We 
used this simplified procedure because a direct integration of the 
model over the stellar surface is a very time-consuming process 
under our approaches. The most general and precise way to pro- 
duce a disk-integrated model spectrum is to 

i calculate (for given c„ and the polar strength Bp"') a num- 
ber of local model atmospheres and synthetic spectra for all 
surface areas from the magnetic equator to the pole, 

ii integrate the local spectra over the visible stellar hemisphere 
for a set of rotational phases and given orientation of the stel- 
lar rotational and magnetic axis. 

At the present stage of our studies we are not yet ready to carry 
out such computations. Thus, our basic modelling strategy is to 
replace the flux from the disk-integrated model which should be 
computed for each observed phase with the flux from a single 
atmosphere model from LLmodels, with parameters which rep- 
resent approximately the disk-averaged means. 



On the other hand, our additional study revealed that even the 
strongest possible intensities of the Lorentz force do not change 
significantly the limb darkening in computations with the mag- 
netic model compared to the non-magnetic one. Therefore, for 
the first investigation we used the above simplified approach. 
This approximate method essentially assumes linear response of 
the properties of magnetic atmospheres to changes in the input 
parameters. 

Finally, the Lorentz force in the cases of dipole and 
quadrupole magnetic field geometries can be written as 



d-Pto tal 

dr 



1 



= ~Pg- 



-(A ± )c 



((pi^ + ^wj^ + ag )). (9) 



where braces mean the disk-averaging operation ((A±) here 
means the "effective" conductivity calculated with the magnetic 
field strength averaged over the stellar disc). Equation (O can be 
written in a simplified form by introducing the effective acceler- 
ation geff as a sum of gravitational and magnetic accelerations so 
that 



(10) 



Results of the disk-averaging procedure of the magnetic pa- 
rameters are illustrated in Table [3] where we adopted fi q /B<j = 
0.0 (dipole) and -2.0 (dipole+quadrupole). Several sets of mod- 
els with the outward and inward directed Lorentz force and 
B q /B ( j ranging from 0.5 to -2.0 with a step of 0.5 were then 
calculated. The best fit to the shape of the observed hydrogen 
line variation was finally found for the following parameters: 
c 2 /ci = 2.5 and B q /B d = -2.0. 

Numerical tests showed that, in order to reproduce the am- 
plitudes of the observed standard deviations due to the profile 
variations in case of the inward directed Lorentz force, the ef- 
fective electric field should be c\ — 2 x 10~ 10 CGS units for 
the dipolar configuration and c\ = 1 X 10~ 10 CGS units for 
the dipole+quadrupole model. In case of the outward directed 
Lorentz force these values should be c\ = 3.5 x 10" CGS 
units and c\ - 1 x ll" 11 CGS units, respectively. Increasing 
c\ further leads to unstable solution. This occurs for both in- 
ward and outward directed Lorentz forces. In the latter case this 
means that magnetic force directed outwards becomes stronger 
than the gravitational force, numerically forcing P gas < 0. For 
example, the critical value of the effective electric field is c\ as 
2 x 10" CGS units for dipole+quadrupole configuration. In 
spite of the fact that inward directed Lorentz force increases 
pressure of stellar plasma, the enormous increase of c\ causes a 
failure of some numerical algorithms implemented in the model 
atmosphere code (mainly hydrostatic equation). This also causes 
numerical problems with the interpolation of partition functions 
for some elements. The estimated critical value for inward di- 
rected force is ci > 2 x 10~ 10 CGS units for dipole+quadrupole 
configuration. 

Figure [7] illustrates resulting dependence of on the 
Rosseland optical depth in the atmosphere of 8 Aur (B q /Z?d = 
-2, c<ilc\ = 2.5). The formation region of the central parts of 
Balmer lines (logTR OSS < -1 and above) demonstrates strong 
changes in the effective gravity during rotation that results in the 
Balmer line variations. Examination of the purely dipolar model 
gives similar results. 

As follows from Eq. ©, the resulting disc-averaged atmo- 
spheric perturbations due to the Lorentz force depend on the 
competition between the effective electric conductivity (A±) and 
the product ^P l „ ■ 2 &g) averaged over the disc. The elec- 
tric conductivity (A ± ) across the magnetic field lines of force 
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Table 3. Results of the surface averaging procedure for the 
dipole and dipole+quadrupole (B q /Z?d = -2.0, Cilc\ = 2.5) 
magnetic field configurations. Field strengths are given in Gauss. 



Rotation 
Phase 


dipole 


dipole+quadrupole 


(B) 


(P[Bg) 


<Z s ( '"> 




0.056 


856 


436 


2158 


-2652 


0.061 


854 


438 


2147 


-2636 


0.091 


842 


449 


2077 


-2533 


0.111 


833 


458 


2016 


-2449 


0.121 


829 


463 


1982 


-2406 


0.250 


812 


480 


1530 


-1913 


0.331 


863 


429 


1332 


-1836 


0.356 


886 


406 


1289 


-1878 


0.377 


905 


388 


1261 


-1931 


0.393 


920 


373 


1240 


-1987 


0.456 


965 


331 


1188 


-2198 


0.496 


975 


321 


1177 


-2256 


0.562 


955 


340 


1199 


-2145 


0.681 


852 


439 


1358 


-1827 


0.775 


806 


485 


1613 


-1989 


0.837 


813 


479 


1833 


-2224 


0.861 


821 


470 


1921 


-2329 


0.936 


853 


439 


2142 


-2627 


0.950 


858 


434 


2169 


-2668 


0.995 


865 


427 


2210 


-2734 
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Fig. 7. Run of the effective acceleration with Rosseland optical 
depth in the atmosphere of 8 Aur obtained for all 20 phases of 
the star's rotation. The calculations were carried out for inward 
directed (upper plot) and outward directed (lower plot) Lorentz 
force for dipole+quadrupole magnetic field configuration. 



is reduced in comparison to the non-magnetic conductivity by 
the full disc-averaged magn etic field (Y, B_y) (the ma g netore - 
sistivity effect, see § 2 in iPikelnerl dl966l) : ICowlingl dl957l) : 
Schluterl (QUI). In our dipolar magnetic field model (B) (and 
the magnetoresistivity as a result) varies with an amplitude 
of about 10% about the average value, which is almost twice 
smaller than the variation of the weighted-average of the hor- 
izontal magnetic field (see Table [3j. Therefore, in case of the 
dipolar geometry, variation of the Balmer lines is produced 
mainly by phase variation of the horizontal magnetic field. In the 



dipole+quadrupole configuration however, the variation of con- 
ductivity also plays noticeable role due to significant rotational 
modulation of (2 B (n) ) (Table [3]), which explains the influence 
of magnetoresistivity on the Balmer line variations for this ge- 
ometry. 

4.5. Comparison with the observations 

We compared observations and theoretical predictions for each 
of the rotation phases. To minimize effects of systematic errors in 
the modelling procedure, we have examined residual theoretical 
and observed Balmer lines, which are obtained by subtracting a 
spectrum at a reference phase {<p = 0.056 where the Balmer pro- 
files have the largest widths) from all the other spectra. Figures[8] 
|9]and[T0]illustrate residual Ha, HB and Hy line profiles for each 
of the observed rotation phases. The positive sign of the residu- 
als implies that the lines at the current phase are narrower than 
those obtained at the reference phase. Here we present only those 
calculations, which are able to reproduce the observed positive 
deviations with satisfactory accuracy. This is achieved for the 
purely dipolar geometry of the magnetic field with the inward 
directed Lorentz force (dotted lines in the figures) and by ex- 
amination of the dipole+quadrupole geometry with the outward 
directed force (thick solid lines). 

The most noticeable effect is seen at phases between 
<p = 0.250 and <p = 0.6. However, both models fit the data 
only in the phase region between (f> = 0.45 and <p = 0.6. We at- 
tribute this disagreement to our assumption of the axisymmetric 
magnetic field geometry (in reality the field can be decentered or 
distorted). Despite this discrepancy we suppose that our model 
reproduces general phase dependence found in observations. 

Comparing results of the fit obtained from both models we 
conclude that observations are better reproduced in the case of 
the outward directed Lorentz force for the dipole+quadrupole 
magnetic field geometry. As can be seen from Fig. [9] [10] the 
inward directed force gives too wide line wings and is unable 
to describe the line cores of the observed residual spectra. In 
contrast, the outward directed Lorentz force shows better agree- 
ment with observations. Such a difference between the two mod- 
els results from the fact that at the upper atmospheric layers of 
the Balmer line formation regions the inward directed Lorentz 
force changes the pressure-temperature balance more effectively 
than the outward directed one (see Fig. |7J. Increasing g e ff, the 
inward directed Lorentz force makes the Balmer line profiles 
non-realistically wider than the observed ones. This makes the 
outward directed magnetic force more preferable in our analy- 
sis. 

The best fit obtained with B q /B& = -2.0 and Ci/c\ = 2.5 
agrees with the longitudinal field modelling results. The out- 
come of the Hy line profile analysis (Fig. [ToT > is in satisfactory 
agreement with the H/3 line modelling. In the case of the Ha line 
(Fig. [8} it is difficult to distinguish predictions of the two models. 

Now, before general discussion, we would like to clar- 
ify some points in the above modelling approaches. The ob- 
served changes in Stark profiles of the Balmer lines demon- 
strate monotonic single-wave variation during full rotation cy- 
cle of 6 Aur. Originally, however, we expected another picture 
of the variations due to predicted dominant dipolar component 
of th e magnetic field in this star. According to previous stud- 
ies dValvavin et al.l 120041 and references therein) we expected 
that if a magnetic star with centered dipolar field shows the ob- 
server both positive and negative parts of its magnetosphere dur- 
ing rotation (that takes place in 6 Aur), the hydrogen lines should 
demonstrate double-wave variations during full rotation cycle. 
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Fig. 8. Residual profiles of the Ho- line for each of the ob- 
served phase relative to the phase 0.056. Thin solid lines are 
observations, dotted lines are theoretical profiles calculated 
for the dipole magnetic field configuration (inward directed 
Lorentz force), thick solid lines are profiles calculated for the 
dipole +quadrupole magnetic field configuration with B q /Ba = 
-2.0 and C2/c\ = 2.5 (outward directed Lorentz force). The 
residual spectra for consecutive phases are shifted in the vertical 
direction. The thin dashed line gives zero level for each spec- 
trum. 



To resolve this difficulty, we have included the quadrupole com- 
ponent, which is still not well-established for 9 Aur, but is plausi- 
ble for the reasons mentioned above. In the contrast to the purely 
dipolar magnetic field, the combination of dipole+quadrupole 
enables us to derive a satisfactory fit to the observed behavior 
of the Balmer line variations. As follows from Table [3j we are 
able to obtain a configuration of the dipole+quadrupole mag- 
netic field for which the maximum visible horizontal component 
is located at those phase points where the varying longitudinal 
magnetic field exhibits extremum (not at the crossover). Such 
a configuration gives the single-wave variation of the Lorentz 
force and reproduces our observations. 

Another important point is related to different orientations 
of the Lorentz force obtained for the two examined geometries 
of the magnetic field. For both geometries with the opposite- 
oriented Lorentz force the sign of the resulting Balmer line resid- 
ual variations is unchanged at the extremal point of the phase 
curve (4> ~ 0.5). In this connection we note again (see above), 
that the Stark widths variations in the purely dipolar model are 
mainly produced by the horizontal magnetic field that provides 



i i I i i i i I i i i i I i i i i I i 




-0.6 



-0.7 



-0.8 



I I I I 

0.061 



0.1 1 1 

0.121 



ffW^At o 




0.331 



^WrU ' V 0-376 



0.393 
456- 
0.496 
0.562 
0.682 




rtty^.^r^fr * ^^ ' ^ 0,775 



0-837- 




t h^Mp A k 0.861 
0.936 



^ *V* H ^ ^M*v > 0.995- 



-30 



-20 



-10 



10 
AA. (A) 



20 



30 



40 



Fig. 9. Same as Fig. [8]but for HB line. 



the minimum inward directed Lorentz force at this phase due to 
the minimum of visible horizontal component (see Tabled. The 
Stark widths of the lines are reduced at this phase in comparison 
to the other phases. In contrast, the dipole+quadrupole geometry 
of the magnetic field provides significant horizontal component 
at this point and minimum magnetoresistivity that results in the 
maximum outward directed magnetic force term. As a result, the 
Stark widths of the Balmer lines are also reduced at this phase 
relatively to all the other phases of the stellar rotation even im- 
plementing the opposite orientation of the Lorentz force. 

Despite the fact that a fairly good agreement with obser- 
vations is obtained, we do not claim that our geometric fit is 
absolutely correct. For example, our assumption about essen- 
tially poloidal centered magnetic field may be wrong. Intuitively 
we may assume, that decentered dipole may also produce the 
Balmer line profile variations which can be similar to the 
observed ones. Besides, we do not exclude the presence of 
more complicated dynamical, non-evolutionary processes such 
as meridional circulation which can also produce similar ob- 
servable features. The problem is too complicated, and the an- 
swer about the true model parameters (also about the true na- 
ture of the Lorentz force) can be obtained only when the ex- 
act geometry of the surface field in 8 Aur will be established. 
In this study we only demonstrate that the observations can be 
described with simple geometrical approaches about magnetic 
field and induced electric currents. We therefore conclude that 
the main achievement of this paper is not our explanation of the 
behavior of the phase-resolved Balmer line profiles, but our ar- 
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Fig. 10. Same as Fig.[8]but for Hy line. 



gumentations for the presence of the detectable Lorentz force in 
the atmosphere of 9 Aur. 

5. Conclusions 

With the aim to probe the presence of the Lorentz force in at- 
mospheres of magnetic CP stars, we have acquired new high- 
precision spectral observations for one of the bright magnetic 
stars, 6 Aur. The hydrogen line profile variations are investigated 
using high-resolution spectra covering the full rotation cycle of 
the star. We have detected and studied variability of the Stark- 
broadened profiles of Ha, H/?, and Hy. Several evidences for the 
presence of significant Lorentz force in atmosphere of 9 Aur are 
found: 

- The characteristic shape of the variation during full rota- 
tion cycle o f the star correspon ds to those describe d by 
iKrollI (119891) and other authors dValvavin et al.1 12004 and 
references therein) as a result of an impact of a substantial 
Lorentz force. 

- Numerical calculations of the model atmospheres with in- 
dividual abundances demonstrate that the surface chemical 
spots can not produce the observed variability of the hydro- 
gen line profiles of 9 Aur. 

- Our model shows good agreement with the observations if 
the outward directed magnetic force is applied assuming the 
dipole+quadrupole magnetic field configuration with the in- 
duced effective equatorial electric field of c\ = lxl0~ n CGS 
units. 



- The dipole+quadrupole magnetic model with the quadrupo- 
lar strength twice stronger than the dipolar one reproduces 
behavior of the phase-resolved H/? and Hy spectra. 



6. Discussion 

In the above considerations we did not examine any mechanisms 
of the magnetic force generation. We have restricted ourselves to 
the purely phenomenological and geometrical description of the 
problem. Let us finally discuss physical processes which could 
create the Lorentz force in the atmospheres of magnetic CP stars. 



6. 1. Secular evolution of the global magnetic field 

It is well known that any variation of the global stellar magnetic 
field (related, for instance, to a global field decay) leads to the 
development of the induced electric currents in conductive at- 
mospheric layers. The Lorentz force, which appears as a result 
of the interaction between the magnetic field and electric cur- 
rents, may affect the atmospheric structure and influence the for- 
mation of spectral lines causing their rotational variability. The 
sign of the observed Lorentz force (outward or inward directed) 
is also found to be important. For example, following IWrubell 
(fl952h . who was the first to present numerical calculations of 
the decay of poloidal interior magnetic field, the decay-induced 
electric field E& have essentially azimuthal form: 



£ ~ g(r)P l M, 



(11) 



where g(r) is a scalar radial function of a distance r from stel- 
lar center (r = [0,/?*]). Considering only the dipole component, 
the induced currents do not change sign at the stellar surface, 
achieving maximum strength at the magnetic equator (where the 
magnetic field lines of force is horizontal) and vanishing at the 
magnetic poles (9 = 0°, 180°). In this case the Lorentz force cor- 
responds to the case of the outward directed decay-induced force 
term at any non-polar area of the stellar surface. Therefore, deter- 
mination of the sign of the Lorentz force provides an important 
constraint which is able to restrict di rection of the field evolution. 

According to Landstreet (1987), who also considered a de- 
cay of an essentially dipolar fossil magnetic field, we should not 
expect any signatures of the Lorenz force in atmospheres of CP 
stars. In his model magnetic field has dipolar geometry through- 
out a magnetic star and electric field induced at the equator can 
be approximated with the following expression: 



eq 



Ct 



(12) 



where is the stellar radius, Z? eq is the strength of the dipole- 
like surface field at the magnetic equator and t is the character- 
istic decay time of the magnetic field. For a typical AO mag- 
netic star with the surface field of 10 kG this leads to the value 
£eq ~ 10~ 13 CGS, which is about 2-3 orders of magnitude 
smaller compared to the strength of the electric currents tested in 
our model calculations and inferred from the variations of hydro- 
gen lines. The former theoretical estimate of E ecj is based on the 
10 10 yr characteristic decay time of nearly dipolar fields. Thus, if 
the observed strong electric currents in the atmosphere of 9 Aur 
have the evolutionary nature, we may conclude that the general 
concept of a slow decay of the fossil, essentially dipolar, fields 
requires for revision. One should examine more realistic struc- 
tures of the internal fields, and this may lead to faster field evolu - 
tion. As it was recently shown by Braithwa ite & Spruitl (|2004), 
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dynamical stability of the global stellar magnetic fields is en- 
sured if the interior field configuration includes both poloidal 
and toroidal components. Diffusive evolution of such a twisted 
field structure may lead to relatively rapid change of the surface 
field intensity. This, in turn, induces a noticeable atmospheric 
Lorentz force during certain periods in the star's life. In this con- 
text advanced evolutionary stage of 8 Aur is remarkable. 

Thus, large amplitude of the induced electric currents may 
indicate that the poloidal magnetic geometry, dominating at the 
surface, becomes significantly distorted inside a magnetic star. 
Such a distortion is very likely to be related to dynamical pro- 
cesses of the field evolution and, possibly, interaction with the 
core convective zone or differential rotation. 

6.2. Non-evolutionary surface effects 

Generally, any global magnetic topology gives a large collec- 
tion of possibilities of generation of the surface electric cur- 
rents even without invoking the global magnetic field evolution. 
There are mechanisms which may produce significant atmo- 
spheric currents even in constant magnetic fields. For instance, 
Peterson & Theysl (1 198 lh considered interaction of the horizon- 
tal component of the magnetic field with a flow of charged par- 
ticles, drifting in the atmosphere under the influence of the radi- 
ation pressure. The result of such an interaction is that drifting 
particles acquire some horizontal velocity component. This leads 
to appearance of the Lorentz force which may be significant for 
hot sta rs (r eff > 18 000-20 000 K). More recentlv lLeBlanc et alj 
(1994) studied a similar physical situation in the context of am- 
bipolar diffusion. Furthermore, interaction of the magnetic field 
and stellar rotation may induce dynamical processes leading to 
the development of additional plasma flows (such as meridional 
circulations) which are able to create electric currents and the 
Lorentz force. 

We have obtained our results under the assumption that the 
equatorial induced electric field is constant in vertical direction. 
This should be close to reality if we deal with evolving global 
magnetic field, whose variation is small throughout the observ- 
able atmospheric layers. This assumption, however, may not be 
applicable if the Lorentz force is related to, for example, the am- 
bipolar diffusion. Nevertheless, our results provide an observa- 
tional test to distinguish these effects. 

As it was shown bv lLeBlanc et alj d 19941) . the ambipolar dif- 
fusion is small for the weak field stars and increases in stars 
with the surface magnetic fields of about 10 kG and higher. 
It should be noted, however, that we have considered a rel- 
atively weak-field star (B « 1 kG). Moreover, the inward di- 
rected Lorentz force predicted by the ambipolar diffusion the- 
ory is not supported by our observations, as we have found a 
much better agreement for the outward directed magnetic force. 
Nevertheless, we do not claim that our approach is the only way 
to address the problem of the generation of the Lorentz force 
in the atmospheres of CP stars. Before making final conclusion 
about the nature of the Lorentz force found in 6 Aur, alternative 
models should be applied to interpret our observational findings. 
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